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A  STUDY  OF  THE  MAXIMUM  SOLID  ROCKET  ENGINE  THRUST  NOZZLE 
SURFACE  FOR  ONE- DIMENSIONAL  TWO-PHASE  FLOWS, 


by  Luo  Junyang  and  Sun  Yue 


Abstract 


This  paper  uses  the  Pontryagin  Maximum  Principle  to  study 
the  problem  of  the  optimal  nozzle  surface  for  one-dimensional 
two-phase  flows.  Choosing  the  slope  as  the  control  variable, 
we  verified  that  the  optimal  surface  was  a  conical  surface 
composed  of  a  group  of  straight  lines  with  the  slope  being 
constant.  Furthermore,  we  discussed  the  method  of  defining  the 
optimal  contour  in  the  subsonic  and  supersonic  regions.  This 
paper  can  provide  a  theoretical  basis  and  computation  reference 
for  the  design  of  nozzles. 


Basic  Symbols 

p  =  pressure 

u  =  speed 

T  -  temperature 

G  =  weight  flow  volume 

=  gas-specific  heat  ratio 

CP9,Cpl  =  ^as'  Part^c^-e  specific  heat 

g  =  gravitational  acceleration 

x  =  x  direction  coordinate 

y  =  radius  of  corresponding  independent  variable  x 

Mp  =  a^2°3  roaterial  density 

s  gas  constant 

=  mechanical  equivalent  of  heat 

=  heat  transfer  correction  coefficient 
=  drag  correction  coefficient 
=  gas  kinematic  viscosity  coefficient 
*  Prandtl  number 
=  particle  radius 


*  interval  of  convergence  half  angle 
P  =  diffusion  region  half  angle 

p  =  external  atmospheric  pressure 

cl 

Lower  Symbols 

g  =  gas  phase 

p  =  particle  phase 

I.  Preface 

Early,  in  1963,  F.E.  Marble  [1]  used  the  one-dimensional 

flow  minimum  hysteresis  method  for  the  design  of  the  optimal 

nozzle  surface  for  two-phase  flows.  He  only  attained  an  integral 

expression  but  this  method  is  not  appropriate  in  engineering 

technology.  In  1968,  A.N .  Krayko,  V.K.  Starkov  and  L.Ye. 

•  *  • 

Sternin  [2]  used  the  supposition  of  the  one-dimensional  flow 
and  employed  the  calculus  of  variations  to  compute  seven  maximum 
thrust  nozzles  in  the  supersonic  region  by  means  of  an  electronic 
computer.  Its  surface  is  a  linear  region  yet  this  method  cannot 
use  an  analytical  solution  to  verify  the  optimal  nozzle  surface. 
From  1966  to  1972,  Hoffman  et  al  [3]  used  the  two-dimensional 
flow  theory  and  the  calculus  of  variations  to  compute  the 
optimal  nozzle  surface  in  the  supersonic  region.  Because  at  pre¬ 
sent  there  is  still  no  perfect  method  for  the  flow  of  the  two- 
dimensional  subsonic  and  supersonic  regions  in  two-phase  flows, 
difficulties  occur  even  if  this  method  is  practically  applied. 

This  paper  investigates  the  problem  of  the  optimal  rocket 
engine  surface  for  one-dimensional  two-phase  flows  using  the 
maximum  principle  to  find  the  best  system  of  solutions.  If  the 
boundary  conditions  are  given,  we  can  find  the  corresponding 
optimal  nozzle  contour  lines. 

II.  Basic  Theories  and  Hypotheses 

1.  Based  on  reference  materials  [4] ,  there  are  the 


following  hypotheses: 

a.  The  flow  is  one-dimensional,  and  stable; 

b.  The  system  exchanges  with  the  external  non- thermal 
quantity  and  the  mass  ; 

c.  Aside  from  the  drag  effect  produced  on  the  particle,  the 
gas  is  considered  non-viscous; 

d.  The  energy  exchange  between  the  gas  and  particles  is 
only  limited  to  the  convection  heat  transfer; 

e.  The  gas  flow  is  the  condensed  flow  of  an  ideal  gas  and 

the  particles  take  on  a  liquid  state  (if  below 

the  solid  point,  it  also  assumes  a  supercooling  liquid 
state) ; 

f.  The  shape  of  the  particles  is  circular; 

g.  The  internal  temperature  of  the  particles  is  uniform; 

h.  The  Brownian  motion  does  not  affect  the  pressure; 

i.  The  area  occupied  by  the  particles  can  be  neglected; 

j.  There  is  no  collision  between  the  particles  and  the 
diameters  of  the  particles  do  not  change  during  the 
flow  process. 

2.  Basic  Theories 

Based  on  the  maximum  principle  of  P.S.  Pontryagin  [5],  we 
used  the  Lagrange  optimum  performance  index  functional  form  of 

/  «  f  1  ♦  («t,  w,  x)dx  (  1  ) 

J  ** 


(1)  Equation  of  state 

*iL  •  „  x 

TjT  u -/(«,  u>,  x)  (2) 

In  the  equation:  ur-(«„«„ . ,«.)  is  the  state  variable  and 

. ,to,)  is  the  control  variable. 

(2)  Conjugate  equation 

dH  ,  \  _  LL rCu(*),  w(»),  *) .  ,  v  *4>C «<(*).  JM  , 

du  A  iu  AV.*;+  -  (J) 
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(3)  H  function 


to(*),  x)+  Xr  (*)/  [«(*),  to(*),  *)  (4) 


(4)  Optimum  principle 

If  to  is  one  inner  point  of  'll)  ,  then  we  have 


iH 

dW 


0 


(5) 


3H 


If  co  Si CO  is  an  arbitrary  point,  then  possibly  does  not 
exist.  If  we  let  H  be  the  minimum  (maximum)  control  Co  ,  and  the 
arbitrary  adjacent  control  is  A>  ,  then 


H(u,  to,  X,  *)-  H(u,  to,  X,  x) 
«6t o 


or 


H(u,  to,  X,  (u,  to,  X,*)  (6) 

(5)  Transverse  conditions 

The  transverse  conditions  are  not  sufficient  to  compensate 
for  the  boundary  conditions.  The  end  point  free  transverse  con¬ 
ditions  are 


»<*<>-(•&),  +(-£sr),* 

N  Cu(*, 0 

Hiu  (*/),  to(*;),  X  (»/),*,]  +  +  *2--v  *0 


In  the  equation:  N  [u(x^) ,xf]  is  the  flow  equation  and  v  is 
the  multiplier. 

Equation  (7)  is  a  two  point  limit  value  problem  which 
requires  the  use  of  the  iteration  method  to  find  its  numerical 
solution.  This  is  naturally  very  difficult.  Yet,  for  one-dimen¬ 
sional  two-phase  flow  engines,  the  equation  of  state  and  value 
functional  both  have  specific  forms,  that  is,  contain  the  linear 
item  of  slope  y.  This  allows  us  to  possibly  use  the  Pontryagin 
Maximum  Principle  to  find  the  theoretically  optimal  nozzle 
surface. 

III.  Equations  of  State 

Based  on  the  above  assumptions,  we  can  obtain: 

1.  The  gas  motion  equations 

(«) 

-1_  dp,  kgdPt  _  APpB 

Pm  dx  Pgdx  pjt,  t 9 ' 
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2.  Particle  motion  equations 


AC  (T,-Tt) 


Gp-gPwy1 


After  equations  (8)  to  (13)  undergo  mathematical  processing 
they  can  finally  change  into 

Letting 


_ zkg pi  • 


®, *» ^4— — f —  (tl,  —  B#) 

g*y 

-~r 


We  finally  obtain 


/ 


dut 


(15) 


From  equations  (]5)  and  (8)  we  can  solve 


(16) 


From  equation  (12)  it  can  change  into 


dT„  2  ACv, 

dx  3  CPfu, 


(17) 


When  equations  (11)  ,  (15)  ,  (16)  and  (17)  are  equations  of 
state  for  two-phase  flows,  they  can  be  written  as  general  forms 


•^--/i(«f,  Pg,u„ 7*„ y,y)-~fl  ( u,y ) 
-^— /*(« t>  />«.  “p,  T„  y,  y) -/,  ( u ,  y) 

■f5f-/.(«»*)-/.  («) 

P**  «p.  p,  y)  “ACu  ) 


We  can  see  from  the  four  equations  of  state  above  that 

T 

aside  from  containing  five  state  variables,  u  =(u  ,p  ,u  ,T  ,y) , 

g  g  P  P 

they  also  contain  a  derivative  of  a  state  variable.  In 
y=g^,  y  is  selected  as  the  control  variable.  Therefore,  there 
are  five  equations  of  state  altogether. 

IV.  Value  Functional 

Because  the  nozzle  is  composed  of  the  convergence  and  expan¬ 
sion  regions  as  shown  in  fig.  1,  the  value  functional  is 
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*1 


,***»«■  ( 4 ) 


(3) 


Fig.  1  Diagram  of  Nozzle  Surface 


Key:  1.  Nozzle  inlet  section 

2.  Nozzle  convergence  region 

3.  Nozzle  expansion  region 

4.  Nozzle  throat  section 

5.  Nozzle  outlet  section 


V.  The  Hamilton  Function 

Two  end  points  are  given  wherein  the  terminal  end  point  can 
be  stable  as  well  as  free.  When  there  are  no  other  constraining 
conditions,  the  H  function  can  be  written  as 


2 

-^^-  +  X>  +  2*( pt-p.)yy 


(19) 


Equation  (21)  is  written  as  the  general  form 


0t(u*X)  +  At(u,  X,y) 


d\„ 
dx  ~ 


AGt 

gny 


I 


ACvt 

Cp,»} 


Equation  (22)  is  written  as  the  general  form 


d\. 


dx 


‘St (u,  X) 


AC 


Equation  (23)  is  written  as  the  general  form 


0*(u,  X) 


2<%-£Ajf.+J?Au. 

gny*  y 


In  the  Equation 


Equation  (24)  is  written  as  the  general  form 

,mS>t(*t,  X)  +  ht(u,  X,  y) 

Vll.  Defining  the  Optimal  Contour  Line 

We  obtain  from  the  equation  of  state  and  proposed  value 
functional  equation  (18)  that  the  H  function  is  the  linear 
function  of  y 

H—H^u,  X)+Ht(u,  X)y 

If  we  first  consider  the  subsonic  region,  our  given  limit 
is 

(25) 

The  aim  of  this  article  is  the  design  of  an  optimal  surface 
and  therefore  y  acts  as  the  control  variable.  Based  on  the  max¬ 
imum  principle,  we  should  select  control  variable  y  so  that  H 
function  takes  on  a  minimum  value.  Because  of  this,  control 
variable  y  should  take  a  numerically  allowable  minimum  value  and 
the  symbol  should  be  opposite  that  of  Hj  (Hj^O  is  equal  to  the 
subsonic  region,  0  is  equal  to  the  supersonic  region) .  Then 


In  this  way,  we  obtained  the  optimal  solution  group:  i.e. 
limited  the  group  of  straight  lines  of  the  slope.  It  is  very 
noticeable  that  the  obtained  maximal  solution  group  can  satisfy 
the  boundary  conditions  of  the  two  points.  From  equation  (26) 
we  can  obtain 


y^-^x +ct 


and  c 2  are  both  undetermined  constants  which  can  be 
determined  by  the  two  boundary  conditions.  For  example,  by  means 
of  the  two  fixed  points  M0(xQ,yo)  and  M^(x^,y^)  or  by  the 
MQ(xo,yo)  point  and  the  convergence  region's  half  angle  A. 

The  first  situation  is 

yi-y^-k  ,(*,-*,) 

ct  -  y» + *.*•  ’•y,  — *• 

xi  x% 

As  a  result,  we  fix  constants  and  c^. 

The  second  situation  is 

tg<* 


Then  we  obtain 


*i-tga 

yo“~*ttga  +  Ci 

**-*+*, tg  a 


That  is,  the  given  boundary  conditions  can  determine  a 
subsonic  region  maximal  contour  line. 

In  the  same  way,  control  variable  y  is  limited  by  the 
following  in  the  supersonic  region 

kt>y>& 

Then,  from  the  maximum  principle  we  obtain 

y-kt 

The  maximal  solution  group  is  also  the  slope  sustaining  the 
limited  group  of  straight  lines.  By  integration  we  obtain 

J-kiX  +  Ct 

It  can  also  suit  the  boundary  conditions  of  the  two  points. 
As  a  result,  the  obtained  optimal  supersonic  region  surface  is 
also  a  ;  >up  of  straight  lines.  In  this  way,  we  can  obtain  the 
optimal  nozzle  surface  as  shown  in  fig.  2. 


Fig.  2  Chart  of  the  Maximal  Nozzle  Surface  When  the  Given 
Boundary  Conditions  are  yQ<,  and 

In  the  aforementioned  value  function 
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r*'  rxt 

/"  J  2*0>,-.p.)<yy</jf+ J  2n(pl-pm)yydx  (18) 

f*‘ 

the  first  item  J  2n(P*~ PJyM*  is  the  drag.  Based  on  the  max- 

*• 

imum  principle,  the  y  selected  in  this  area  is  the  minimum  value 
ai^i  thus  the  entire  drag  is  minimum;  the  latter  item 
j*  2*0*,— Pm)  is  the  thrust.  Based  on  the  maximum  principle, 

the  y  taken  in  this  region  is  the  maximum  value.  Thus ,. maximum 
thrust  is  produced  in  this  area  and  from  the  sum  of  the  two 
items  we  then  obtain  the  maximum  thrust  surface.  The  reason  we 
can  obtain  this  result  is  because  the  Hamiltonian  function  H  is 
the  J  linear  function.  Based  on  the  Pontryagin  Maximum  Principle, 
after  y  takes  the  maximum  or  the  minimum,  the  value  function  is 
also  the  integration  of  the  thrust  and  drag  and  takes  the 
maximum  value.  This  is  the  physical  significance  of  this  paper 
concerning  the  optimal  surface. 

Vlll.  Conclusion 

The  maximum  thrust  nozzle  surface  for  one-dimensional  two- 
phase  flows  is  a  group  of  conical  surfaces  when  the  slope  is 
limited.  If  we  give  the  boundary  conditions,  the  optimal  nozzle 
surface  is  a  conical  surface.  Thus,  we  obtain  a  deduction;  when 
designing  a  one-dimensional  two-phase  flow  nozzle  surface,  con¬ 
vergence  half  angle  of  the  subsonic  region  should  be  made  as 
small  as  possible  and  expansion  half  angle  of  the  supersonic 
region  should  be  made  as  large  as  possible. 
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ITERATIVE  CALCULATION  METHOD  FOR  THE  EQUILIBRIUM  COMPOSITION 
OF  THE  COMBUSTION  PRODUCTS  OF  COMPOSITE  PROPELLANTS* 

by  Wan  Junhua 


Abstract 


This  paper  presents  a  simplified  calculation  method  which 
solves  the  equilibrium  composition  of  the  combustion  products 
of  composite  propellants.  Using  element  mass  conservation 
equations  and  equilibrium  equations  for  chemical  reactions, 
this  paper  gives  the  iterative  calculation  formulae  for  cal¬ 
culating  the  partial  pressure  of  the  gas  composition.  The 
calculation  procedure  is  simple,  clear  and  much  easier  for 
calculating  than  the  concentration  method.  The  calculation 
results  are  very  close  to  those  obtained  by  the  concentration 
method,  the  partial  pressure  convergence  of  the  composition  is 
very  good  and  to  calculate  4  times  is  generally  very  accurate. 
This  proved  that  the  simplified  calculation  method  proposed 
in  this  paper  is  reliable. 


Main  Symbols 


A,B,C,D,E,F , L 

G 

K 


m 

M 


pi 

P£ 

Q,R,S,U,V 

y 

z 

*"$,r 

Cp , Hp  °p'Np 

* 


calculation  factors 

weight  or  weight  fraction  of  element 

partial  pressure  equilibrium  constant 

gram  molecular  weight  of  composition 

weight  of  composition 

gram  molecular  number  of  composition 

partial  pressure  (atm)  of  i  type 

total  pressure  (atm)  of  combustion  products 

component  constants  of  propellants 

weignt  fraction  or  weight  of  combustion  products 

gram  atomic  number  of  element  in  propellant 

constant  factors 

Sp  element  weight  component  of  propellant 
calculation  factors 


♦Received  Dec.  29,  1979 


I.  Preface 

There  are  various  different  methods  to  calculate  the  equili¬ 
brium  composition  of  the  combustion  products  of  propellants,  for 
example,  the  Huff,  White,  Brinkley  and  simplified  methods  etc. 

The  three  aforementioned  methods  consider  the  possibly  produced 
combustion  products  composition.  Because  the  unknown  numbers 
are  numerous  and  calculation  is  complex,  it  is  generally  neces¬ 
sary  to  use  an  electronic  computer  for  completion.  However,  the 
simplified  method  overlooks  the  small  quantity  composition  con¬ 
tained  in  the  combustion  products  so  that  calculations  are 
greatly  simplified  and  common  calculation  tools  can  be  used  for 
completion. 

This  paper  carries  out  investigations  on  calculation  methods 
for  the  equilibrium  compositions  of  combustion  products  of 
propellants  with  9  elements:  C,H,0,N,Cl,S,Al,Cr  and  Cu.  The 
combustion  products  of  these  types  of  propellants  contain  a  gas 
phase  and  condensed  phase.  After  overlooking  the  very  small 
content  of  composition,  there  are  three  types  of  condensed  phase 
compositions:  Al203,  Cr203  and  Cu.  Based  on  the  chemical  reaction 
equations,  we  used  the  mass  conservation  principle  to  directly 
find  their  contents.  After  using  the  separation  method  of 
independently  calculating  the  condensed  phase  composition,  the 
unknown  number  of  the  gas  phase  system  decreases  to  16  types  of 
compositions.  The  control  equations  of  the  partial  pressure  of 
the  major  compositions  are  still  water  gas  reaction  equilibrium 
equations . 

II.  Equilibrium  Composition  Equations 

1.  The  Equilibrium  Composition  of  Combustion  Products 

Based  on  materials  in  the  references,  the  major  combustion 
products  gas  phase  compositions  of  composite  propellants  have 
C02,  CO,  H20,  H2,  N2,  HC1 ,  S02,  H2S,  OH,  02,  NO,  Cl2,  0,  H,  N 
and  Cl  etc.  and  the  condensed  phase  compositions  have 
A1-0-,  Cr,0-  and  Cu  etc.  Under  high  temperatures,  the  A1_0_ 


separates  into  A^O  and  AlO.  However,  because  polysulfide 
propellants  do  not  have  very  high  energy,  the  temperature  of 
the  combustion  products  is  generally  about  3,000  K.  Under  this 
temperature,  the  A^O^  is  generally  not  able  to  separate  into 
A^O  and  AlO. 


Hypothesis:  all  of  the  sulfur  in  the  propellant  participates 
in  the  chemical  reaction  and  produces  S02  and  H2S  gas. 

Because  the  sulfur  content  in  composite  propellants  is  not 
high,  generally  only  about  7%,  the  oxygen  content  reaches  to 
over  40%.  Moreover,  a  good  deal  of  H2  gas  can  be  produced  in  the 
combustion  products.  Sulfur  easily  oxidizes  and  produces  SC^. 
Under  high  temperatures,  the  sulfur  can  also  directly  combine 
with  hydrogen  to  produce  H2S  gas.  Therefore,  this  hypothesis  is 
close  to  the  actual  situation. 

Based  on  the  above  analysis,  for  the  polysulfide  propellant, 
there  are  19  types  of  equilibrium  compositions  of  this  method  of 
calculation:  02,  OH,  NO,  F^S,  Cl2,  H,  Cl,  0,  N,  N2»  S02,  HC1, 
C02,  CO,  H20,  H2,  Al2°3 ,  Cr2°3  and  Cu* 

The  condensed  phase  composition  is  directly  solved  according 
to  the  chemical  reaction  equations  and  mass  conservation  prin¬ 
ciple.  Below  we  explain  by  using  A^O^  as  an  example. 

Hypothesis:  the  aluminum  and  chromium  in  the  propellant  is 
completely  oxidized  according  to  the  chemical  equivalent  ratio. 

2A1+— |— Ot“AltO» 

H.M  41  101.9* 

M  x\  (A^Oj)aI  ^  K\lO% 
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The  weight  of  the  produced  condensed  phase  A1  0  is 

«  J 


1/  _  191.96  i g  . 

"A,*°s  53^6  Mai  k* 


The  weight  of  the  consumed  oxygen  is 


6Ma' 


In  the  same  way,  it  is  easy  to  find  M_  and  (M  )„  . 

1  Cr-0„  o_  Cr 

2S  2 

Therefore,  after  the  combustion  of  1  kg  of  propellant,  the 

produced  condensed  phase  weight  (or  weight  percentage)  Y  and 

s 

gas  phase  composition  weight  Y^  are 

Y,**Mk  ito,+Afcto,+Woi  kg 
kg 


2.  System  of  Equations  for  the  Equilibrium  Composition 
The  element  mass  conservation  equation:  the  weight  of  a 
certain  element  in  1  kg  of  propellant  should  be  equal  to  the 
weight  of  that  element  in  the  combustion  product.  That  is 

Cr”  ~  2  wp, "  (A=°* + pco)  (  1  ) 

1  y. 

2m,/,  (  2^h,o  +  2pHl+ Atci  +  2p*lS +p»+ pOH)  (  2  ) 

16K 

2miPt  ( ^°*  ^C0*  +  Pco  +  />H,0  +  2^50;  +  .PnO  +  PoH  +  Po)  (  3  ) 

Nr»  i^-(2pNt-t-pN0-*-^)  (  4  ) 


35.457, 

1m,pi 


"CpHCl  +  2Pcij  +  Pci) 


I 
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Because  the  total  weight  of  the  gas  phase  composition  is 

Y  k  ,  the  right  sides  of  equations  < 1 ) — ( 6 )  are  all  multiplied 
y  a 

by  factor  Y  . 

9 

The  dissociation  reaction  equilibrium  equations 


K  i -PcoPo,*'1/ Pcot 

(7) 

Kt  “  pHiPot*' */pHtO 

(8) 

PcoPh  jo/  PcOj^H, 

(9) 

Ki  *  PoHpHi*'*/ PH'O 

(10) 

K***  PnoV PtitPoi 

(11) 

K  ■ = Ph1/ pHt 

(12) 

Kt  —  Po  /  Pot 

(13) 

K,-pHl/pHl 

(14) 

PhPcJ Phci 

(15) 

Based  on  the  dissociation  reaction 

h,s^h,+  sg*)  (1) 

Key:  1.  Condensed 

Therefore 

K***  Phi/Phis 

K  It™  Pci*/ Pc\% 

Based  on  the  Dalton  Law 

(18) 

Among  the  18  above  designed  equations,  because  only  two  of 
equations  (7)-,  (8)  and  (9)  are  independent,  therefore  the 
equilibrium  composition  system  of  equations  only  have  17  in¬ 
dependent  equations.  The  unknown  number  is  also  only  17  (the 


(16) 

(17) 
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total  pressure  of  the  combustion  product  and  the  partial 
pressure  of  16  types  of  compositions) .  Based  on  mathematical 
theory,  this  system  of  equations  is  understandable. 

III.  Simplified  Method  of  a  System  of  Equations 
From  equations  (5)/(4) ,  we  obtain 


_  .35. 45  Png  +  +  Pm 

14  2pNt+pNO+PN 


Then 


P*a  +  2fin*  +  Pa  _  14  Cl,., 

.  2/>Nt+A«>+Ai  35.45  NP  —  1/ “const  (19; 


VC=pHCl+l 


(20) 


In  the  formula 


C  mx2pSt+pNO+pli^2pnt+  A  (21) 

A~  P*o+pn  (22  a) 

L—2pcll+pa  (22  b) 


Therefore 


P  hci  —  CU  —  L 

From  equations  (6)/ (4),  we  obtain 


(23; 


PiO.+pHiS  _  14  Sp 

C  32~  "NT'  V  =CODSt 


(24) 
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Letting 


F  “Ph,* 


Therefore 


Pi  ot 


CV-F 


From  equations  (1)7(4),  we  can  obtain 


Pcot  Pco  7 

C  6 


0  “const 


Therefore 


Pcot  +  Pco  —  CQ=  <t> 


Then  we  obtain 


Pco*  <J>  —  Pco, 


From  equations  (3)/(l),  letting 


D  *  2po,  +  Pno + Poh  +  Po 

% 

we  can  obtain 


2pCO;  +  Pco  ~t~  Ph,0  +  2ptOt  4-  D  3 

Pco*  +  Pco  4 


S  “const 


PcOj  +  (Pco,  +  Pco)  +  Ph,o4"  D  +  2Pw>,  “  S  (Pco,  +  Pco) 
PcOj  +  pHiO+Z)*(S  —  1  )  (Pco,  +  Pco)  ~  2Pso, 


i; 


After  (jp-  +PC0)  is  replaced  by  equation  (27)  and  Pso 
2  2 
replaced  by  equation  (25)  ,  the  above  equation  changes  to 

Pc0i  +  pHto+D-=(S  -  1  )CQ—2CV+2F 

Therefore 

J>h,o=  C(Q(S  —  l  )-2F)-jD+2F -/>«>, 

That  is,  we  obtain 


•PhjO=  X  —pcot 
In  the  formula 


(31) 


X=Ca-£)+2F 
®  “Q(S-  1)-2F -const 

From  equations  (2)/ (4),  we  obtain 

2pttlo+2pH+p»a+  B - 14-H2_»/?.const 

C  Nr 


(32) 

(33) 


(34) 


In  the  formula: 


B  —  i>OH  +  i>H  +  2^H,s  — jPoh+/,h  +  2F  (35) 

2  />MfO +2/>hj  +  /’hci+  B  =*  CF 


Then 

2(i>HjO+/>Ht)“('^?—  Phci~  B  (36) 

When  the  pHC1  is  replaced  by  equation  (23),  equation  (36) 
changes  to 

2(/>H,o+/>H1)“CF-a/+  L-B 

P»i—^-CC(R-U)+L-B}-pHlo 


When  the  P  _  is  replaced  by  equations  (31) ,  (32)  and  (33) 

H2° 

pHt’m~T~^C(R~lJ)+L-B')-tCa-D+zF-pcoj 

pH*=C[-r<-R-V)-a]  +  -^-(L-B)  +  D_-2F+pCOt 
That  is,  we  obtain 

pHt=  V+ficot  (37) 

In  the  equation 

'V“c$+-Y-(L-B)+D-2F  (38) 

P  ~U)—  o.  =  const  (39) 

When  equations  (28) ,  (31)  and  (37)  are  substituted  into 
equation  (9) 

^^m‘(<^>~Fcot)(X  —Pcot)/Pcot(V  +Pa>t) 

After  the  above  equation  is  expanded  and  simplified,  we  can 
obtain 


(K>-  1  )/>*co,  +  (<l>  +  X+'P^,)/>col-<I>X“  0 


(40) 


Equation  (40)  is  the  control  equation  for  controlling  the  major 
composition's  partial  pressure  Pqq  •  A  positive  root  of  equa¬ 
tion  (40) ,  that  is  Pco  •  is  separately  substituted  into  equations 
(28)  ,  (31)  and  (37)  and  so  we  find  the  values  of  Pco  0  an<j 
P 

H2 


2.  The  p„  Formula 


The  calculation  factors  $  ,  X  and  ^  in  equation  (40)  all 
contain  factor  C  and  C=2pN  +A.  Therefore,  before  solving 

equation  (40)  we  must  first  find  the  value  of  PN2. 

Based  on  equation  (18) 


Pc**  Pcot+ Pco+ t>HtO+ pHt+ P*t+ Phc\+  piQl+  E  (41) 


In  the  formula: 

E  ^ Pot  +  Poh  +  Pno+  Pnts  +  ■Pci»  +  /,o  +  .Pci  +  />n  +  ,Ph  (42) 

From  equation  (36)  we  obtain 


2(Ai,o+ .Ph1+.Phci)“C^4‘,Phci~  B 

After  the  P  on  the  right  side  of  the  above  equation  is 
HCl 

replaced  by  equation  (23) ,  it  then  changes  to 

/HijoH"  Phx+ Phci |-C(/?  +  C/) ~(B  +  L)  (43) 


The  P  ^  in 
S02 

(P  .+  P  ) 

co  2  eo 

substituted 


equation  (41)  is  replaced  by  equation  (25) ,  the 
is  replaced  by  equation  (27)  and  at  the  same  time 
into  equation  (43) .  After  merging  we  can  obtain 


Pc-Cy+E — B  +  L)+p*t 


In  the  formula 


7  (R  +  U)+  Q  +  V  “const 

« 

Therefore 


(44) 


Pc-(2fiHi+A)y  +E - |-(B  +  I)+/>Nl 
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Then  we  obtain 


(45) 


,  pc+(B  +  L)/2-Ay-E_ 

Phi  “  2  Y  +  1 


3.  Iterative  Calculation  Method  for  Equilibrium  Compositions 
The  above  derived  formulae  are  arranged  into  two  systems  of 
equations 


^:=[i>c  +  4~(5  +  L)~^Y'£  j/(2  Y  -f  1) 

pc c.«  {-  (<U  +  X  +'VK,')  +  t(4>  +  X 
+  4  (K>-  1  )<t>X)°**)/2(A%- 1) 


pco™  C*  Pc Ot 
X—PcOi 

Ph,  “  W  +  Pco% 
P*ot-CV-F 

pHCl  ™  i 


(46) 


The  set  of  equations  in  (46)  are  equations  for  the  partial 
pressure  of  seven  major  types  of  compositions  of  combustion 
products . 

We  can  calculate  the  set  of  equations  for  the  partial  pres¬ 
sure  of  nine  secondary  compositions  from  equations  (7) ,  (10)  and 
(11) - (17)  . 


Poi^iKiPcot/Pco)1 
Po»™  K»pHio/pHt*m* 

PnO™(.  KtPNtPoi)*'* 
Ph  ™  (  ^*^S.pH|)^** 
Po™  (A,^Oi)6"* 
Ph™  (.KiPst)*'* 

Pc 1  ™  K lpHC\/ pH 

Ph%$  “  Ph%/  K  • 

•PcU^AnV-^lt 


(47) 
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The  iterative  calculation  procedure  is  as  follows:  assuming 
the  partial  pressure  of  the  9  types  of  secondary  compositions 
is  zero,  we  can  then  use  the  set  of  equations  in  (46)  to  directly 
find  the  first  degree  approximation  of  the  partial  pressure  of 
7  types  of  major  compositions.  This  value  is  substituted  into  the 
set  of  equations  in  (47)  to  find  the  second  degree  of  approxima¬ 
tion  of  the  partial  pressure  of  9  types  of  secondary  compositions. 
Then,  we  use  (46)  to  calculate  the  partial  pressure  of  the  major 
compositions.  This  iterative  procedure  is  continuously  carried 
out  until  the  calculation  results  converge  to  the  required 
accuracy. 

The  iterative  convergent  conditions  vary  with  the  required 
accuracy.  If  accuracy  requirements  are  high,  the  convergent 
conditions  are 


pru-pr<  5  X  10-4«tm(A>l  atm) 

pru~pr<  5  xi0'*atmQ>/<l  atm) 


(48) 


If  the  accuracy  requirements  for  the  partial  pressure  are  not 
very  high,  the  convergent  conditions  can  appropriately  widen. 

4 .  Check  Formula 

When  the  two  sides  of  equations  (27)  ,  (32)  and  (38)  are 
increased,  after  arrangement  we  obtain 

<l>  +  X  +  V-CC?+-5-Ci? — l-CU  +  -^ - f- 

The  CQ  on  the  right  side  of  the  above  equation  is  replaced  by 
equation  (27) ,  the  CU  is  replaced  by  equation  (23) ,  the  CR  is 
replaced  by  equation  (36)  and  after  merging  and  arrangement,  we 
obtain 


<I>  +  X  +  ¥  *  Pc  - (Pm + Phci  +  Pso,  +  E ) 


(49) 


Before  solving  the  control  equation ,  we  must  use  equation 
(49)  to  check  as  to  whether  coefficients  ,  X  an<*  ^  are 
correct,  the  value  difference  between  the  two  sides  of  equation 
(49)  is  very  small.  Generally,  it  should  be  smaller  than 
1x10  ^atm  (formed  by  the  calculation  accumulation  errors) . 

The  partial  pressure  values  of  PCq  ,  Pco  and  PH2  should  be 
checked  by  equation  (9)  each  time  of  iterative  calculation.  The 
partial  pressure  of  the  entire  composition  should  be  checked  by 
equation  (18) .  The  final  results  of  the  partial  pressure  of 
each  type  of  composition  must  be  checked  by  mass  conservation 
equations  (1)— (6) . 

IV.  Calculation  Examples 

Example:  the  hypothetical  chemical  formula  of  a  certain  known 
trademark  polysulfide  propellant  is 

Cu.«isoHa«.4»loOjs.,sj0N*.j»„C!».ii7»S..tl|)0Alj.tiT7Cr,.i»„Cu0.I,ii 

With  combustion  chamber  pressure  ^  =60  atm  and  combustion  product 
temperature  ^  =2900K  we  find  the  pressure  of  the  equilibrium 
composition  of  the  combustion  product. 

To  compare  the  results  obtained  from  this  method  and  the  con¬ 
centration  method  12] ,  the  very  small  amount  of  composition 
contained  in  the  combustion  product  was  disregarded,  that  is 
we  assumed 

/,CI|m^HjS“Jp0=IpNs“  0 

1.  The  Weight  Component  of  the  Propellant  Element 

Cp-0. 14020,  Hp-0.03848,  Op  =  0. 38328 
Nr“0. 07353,  C1p»0. 18423,  Sp  =  0. 07296 
A1p-0. 08897,  Crp»0. 00826,  Cuf  =  0. 01009 

Check:  JG.«1. 00000 
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The  oxygen  Og  kg  required  for  producing  the  gas  phase 
composition 


Oj-Of-Ox.-Oc-0. 38328-0. 07914-0. 00381 -0.30033  kg 

2.  Weight  of  the  Condensed  Phase  and  Gas  Phase  Composition 
in  the  Combustion  Product 

•  K,-MAIiO,  +  Mc.,o,+-Mc.-0.  16811  +  0.01207  +  0. 01009-0. 19027  kg 

Kf  — 1.0  — Ki— 0.80973  kg 

3.  Propellant's  Component  Constant  and  Constant  Factor 

0-2.22449,  R  —7.32653,  5  -1  60662,  17- 0.98948 
V  —0.43411,  a -0.48120,  0 -2.68733 
Y -6.81661,  2Y +  1  -14.63322 

4.  Calculation  of  Partial  Pressure  of  Gas  Phase  Equilibrium 
Composition 


The  calculation  procedure  is  as  shown  in  table  1 


)** 

- (3)***a 

tj-UJclT)  "  -- — - . 

Pi 

1 

2 

3 

4 

1 

(KiPco^Pco)1 

Pot 

0 

0.0000266 

00000767 

0.0000267 

2 

KiPU  ;0/Pyl 

POH 

0 

0.02422 

0.02407 

0.02408 

3 

IKiP.UPo^’A 

PNO 

0 

0.00112 

0.00112 

0.00112 

4 

- 

Pm 

0 

0.54247 

0.53857 

0.53859 

S 

(Ktpuci/Pu'i 

Pc 

0 

0.13567 

0.13357 

0.13360 

6 

B  * Po»*Pv 

0 

0.56669 

0.56264 

0.56267 

7 

D  •  2 po4  *  Ps o  *  Poh 

0 

0.02540 

0.02525 

0.02526 

3 

E  «  Pot  *  Pott  *  Pa o  “■  Ph  *  Pci 

0 

0.70351 

0.69736 

0.69742 

1*’ 

_  /£Nattac» 

Pi 

i 

2 

i  3 

1 

4 

9 

Ay 

0 

0.00763 

i  0.00763 

0.00763 

10 

( B  +  pcO/2 

1 

0 

0.35118 

0.34811 

'  0.34814 

11 

Pc-E-Ay+iB  +  PcO/2 

1 

1 

60 

59.64004 

59.64312 

1  59.64309 

12 

Cpc-E-Ay+(B+  pci>/23/2Y  + 1 

1  pH, 

4.1003 

4.0757 

|  4.0759 

4.0759 

13 

C  «  2/>n,  +  A 

8.2006 

8.1525 

8.1529 

8.1529 

14 

4>-CC? 

' 

18.2422 

18.1352 

'  18.1361 

18.1361 

15 

CV 

PSOi  i 

3.5600 

3.5390 

3.5392 

3.5392 

16 

Co 

3.9461 

3.9230 

3.9232 

3.9232 

17 

cp 

i 

22.0377 

21.9082 

21.9093 

21.9093 

18 

CU  -Pa 

PH Cl 

8.1143 

7.9312 

7.9337 

7.9337 

19 

X-Ca-D 

3.9461 

3.8976 

3.8980 

3.8979 

20 

-(B-pcO/2 

0 

-0.2155 

-02145 

1  -0.2145 

21 

V -C0+D-(B -Pcl/t 

22.0377 

21.7181 

21.7201 

21.7201 

22 

<t>+X  +  V<*X)  (7) 

(9) 

44.2260 

43.7509 

43.7542 

43.7541 

23 

Pc- (E  ■*  Ph,  +  Pso4  +  PHCl) 

* 

w 

44.2254 

43.7507 

43.7538 

43.7538 

24 

VK,  +  < t>*X 

180.9038 

178.4466 

178.4623 

178.4622 

25 

32726.1849 

31843.1891 

31848.7925  : 

31848.7568 

28 

4</C,-l)<t>X 

i 

1785.8163 

1753.5237 

1753.7896 

1753.7446 

27 

\WK,  +  <t>+X)l+  4  l>OX>°  * 

185.7741  ! 

183.2941  ! 

183.3101 

183.3096 

28 

c@-  @3/12. 404 

Pco,  i 

0.3926 

0.3908 

0.3908 

0.3908 

20 

♦  -  Pco , 

Pc  0  1 

17.8496 

17.7444 

17.7453 

17.7453 

30 

X  -Pco, 

Ph.o  | 

3.5535 

3.5068 

3.5072 

3.5071 

31 

V  +  pco, 

PH.  | 

22.4303 

22.1089 

22.1109 

22.1109 

32 

Pco’pHjo/pco,*PHt<* X)  (8) 

: 

7.2028 

7.2019 

7.2025  j 

7.2023 

33 

Px-2P, 

pc  i 

60.0006 

60.0002  | 

60.0004  j 

60.0033 

Table  1 


Table  1  Calculation  Procedure  for  the  Partial  Pressure  of  the 
Combustion  Product  CPc=60  atm,  TC=2900K) 

Key:  1.  Number 

2.  Calculation  formula 

3.  Number  of  times  of  Iteration 

4.  Number 

5.  Calculation  formula 

6.  Number  of  times  of  iteration 

7 .  Check 

8 .  Check 

9 .  Check 


Mass  conservation  equations  ( 1) — ( 6 )  are  used  to  check  the 
partial  pressure  of  the  fourth  iterative  calculation. 

12 

«-<  m.p. =1256. 9492,  the  checking  results  are  listed  in  table  2 


(  3)  ft  2 


Table  2 


c 

H 

0.14020 

0.03040 

0.14020 

0.03040 

N 

a 

0. 07353 

0.10423 

0.07353 

0.10423 

Key:  1.  Weight  component 

2.  Element 

3.  Propellant 

4 .  Combustion  product 


We  can  see  from  this  that  the  calculations  are  totally 
correct. 


V.  Concluding  Remarks 

1.  The  Congruence  of  the  Partial  Pressure  of  Each  Type  of 
Composition 


The  major  composition 


A  Pmtx  /*Hi04>  PhiO**  =  0  •  0001 

Calculations  showed  that  the  congruence  of  the  partial  pressure 
of  each  type  of  composition  is  very  good  and  completely  agrees 
with  the  congruent  conditions  of  equation  (48) .  The  congruent 
speed  was  also  very  fast  and  the  fourth  iteration  exceeded  the 
congruent  accuracy. 

2.  Comparison  With  the  Concentration  Method 
The  combustion  product  compositions  in  this  paper  were 
indicated  by  pressure.  The  secondary  composition  formulae  are 
those  of  (47)  .  However,  when  the  concentration  method  is  indicated 
by  the  gram  molecular  number,  then 


«oh- 

By  comparing  the  above  two  formulae,  we  can  see  that  when  the 
composition  is  indicated  by  pressure,  none  of  the  calculation 
formulae  have  the  exponential  factor  (P^/n^)  •  Therefore,  the 

formulae  are  simple  and  calculations  are  convenient. 

Based  on  table  3  we  know  that  the  maximum  corresponding  error 
of  the  partial  pressure  is 

f  Aj>(  1  __4£f!£-« 0.0024 -0.24* 

l  pi  J.U  JpNO 


Therefore,  it  is  said  that  the  results  of  the  two  methods  are 
identical. 


(8) 

*  ft 


bpt/pt 


(9£.  1.  ft£ftC2  3ftS*lft«tl-*i  2.  ftftfttt-*ftil!ftjft»tt£ft'?ft,  «£*.  Pi- ptn,/nt**&.B.))t 
3.  *#»»«*  4 


Table  3  Comparison  of  the  Calculation  Results  of  This  Method 
and  the  Concentration  Method  (Pc=70  atm,  TC=2900K) 

Key:  1.  Partial  pressure 

2 .  Method 

3 .  Concentration  method 

4.  This  method 

5.  Partial  pressure 

6 .  Method 

7 .  Concentration  method 

8.  This  method 

9.  Notes:  1.  Propellant  data  calculations 

based  on  reference  [2] ; 

2.  Composition's  gram  molecular 

number  from  calculations  by 

concentration  method: 

p.=p  n./n  conversion  into 
i  c  i  g 

pressure; 

3.  Both  methods  have  the  results 
of  four  times  of  iteration. 


3.  This  method  is  also  suitable  for  other  types  of  propellants 
especially  for  double-base  propellants  . 
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A  COMMENT  ON  THE  PLATEAU  COMBUSTION  MECHANISM  OF  DOUBLE-BASE 
PROPELLANTS 

by  Ma  Xieqi 
Abstract 

This  is  a  comment  using  58  references  on  the  plateau  combus¬ 
tion  mechanism  of  double-base  propellants.  The  experiments, 
evidences  and  view  points  for  each  principle  are  briefly 
presented,  analyzed  and  reviewed.  Views  concerning  future  work 
are  brought  forth. 

I.  Preface 

The  burning  rate  of  plateau  double-base  propellants  in  a 
certain  pressure  range  appears  to  be  unrelated  to  pressure. 

This  is  very  advantageous  for  rocket  engine  work.  For  this 
reason,  a  large  number  of  scientific  research  personnel  have 
attempted  to  clarify  the  plateau  combustion  mechanism  so  as  to 
better  control,  utilize  and  develop  plateau  propellants.  This 
paper  will  briefly  analyze  and  comment  on  the  various  view 
points,  evidences  and  experiments  related  to  the  plateau  com¬ 
bustion  mechanism  of  double-base  propellants. 

II.  The  Subsurface  Condensed  Phase  Theory 

Camp  et  al  [1]  (1958)  pointed  out  that  after  a  lead  compound 
is  added  into  a  double-base  propellant,  the  lead  compound  can 
absorb  the  ultraviolet  rays  emitted  in  the  flame  region.  This 
changes  the  ester  nitrate  degraded  photosensitizer.  As  a 
result,  this  accelerates  the  ester  nitrate  degradation,  raises 
the  burning  rate  of  the  propellant  and  produces  super  burning 
rate  combustion.  Yet,  following  the  rise  of  the  propellant's 


burning  rate,  the  vanishing  speed  of  the  condensed  phase 
quickens,  the  effective  operating  time  of  the  photochemical 
reaction  in  the  subsurface  reaction  region  shortens,  the  level 
of  the  super  burning  rate  combustion  drops  and  plateau  combus¬ 
tion  drops  and  plateau  combustion  is  produced.  When  the 
pressure  increases  further,  the  condensed  group  produced  by 
the  lead  compound  obstructs  the  absorption  by  the  subsurface 
of  ultraviolet  rays.  This  causes  the  burning  rate  to  decrease 
and  produces  mesa  combustion. 

Camp  et  al  [2]  (1975)  used  a  catalyst  with  different 

excitation  modes  radiated  by  ultraviolet  light  to  develop  the 
photochemical  reaction  theory.  Experiments  by  Kubota  et  al  [3] 
proved  that  increases  in  burning  speed  are  related  to  the  amount 
of  light  radiated  from  the  outside  but  is  unrelated  to  the 
amount  of  ultraviolet  light.  However,  Camp  considers  that  the 
photochemical  reaction  theory  is  still  effective  for  producing 
certain  catalysts  with  high  pressure  plateaus.  Under  low 
pressure,  catalytic  combustion  has  an  unluminous  flame  yet  lead 
compounds  can  still  increase  the  propellant's  burning  rate. 
Therefore,  at  least  under  low  pressure,  the  super  burning  rate 
combustion  is  unrelated  to  the  photochemical  theory. 

2.  The  Chemical  Equivalent  Shift  Theory 

Kubota  et  al  [3,4]  (1973)  used  high-speed  photography  and 

fine  thermocouples  to  research  the  reactions  of  combustion  waves 
in  each  region.  Lead  salt  causes  the  combustion  surface  to 
accumulate  carbon  bearing  substances  and  the  corresponding 
burning  rate  also  increases.  However,  the  lead  salt  does  not 
change  the  amount  of  heat  released  by  the  surface.  It  only 
causes  the  hissing  region's  thermal  conduction  feedback  to 
increase  50-100%.  They  consider  that  the  appearance  of  carbon 
is  the  result  of  the  ester  nitrate  decomposition  producing 
changes  and  the  decrease  in  the  produced  quantity  of  aldehyde. 
This  causes  the  ratio  of  aldehyde  and  NC^  to  become  even  closer 
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to  a  chemical  equivalent.  However,  the  reaction  of  aldehyde 
and  NC>2  is  fastest  under  a  chemical  equivalent.  a  result, 
this  increased  the  reaction  heat  and  raised  the  burning  rate. 
Following  the  raising  of  the  pressure,  the  effective  time  of 
the  catalytic  effect  shortened,  the  quantity  of  aldehyde  in-  | 
creased,  the  ratio  of  aldehyde  and  N02  departed  from  the 
chemical  equivalent  and  the  aldehyde  and  N02  reaction  speed 
lowered  which  produced  plateau  combustion.  Kubota  also  set  up  a 
mathematical  model. This  theory  employed  the  research  results  of 
Pollard  et  al  [5]  on  the  propagation  speed  of  the  formaldehyde 
and  NC^  combustion  flame.  When  the  C^O/NC^  ratio  is  1/1.4,  the 
burning  rate  is  maximum.  Its  reaction  is  as  follows: 


This  theory  does  not  concretely  explain  how  the  lead  compound 
changes  the  path  of  ester  nitrate  decomposition.  It  can  only 
explain  super  burning  rate  combustion  and  plateau  combustion  but 
cannot  explain  mesa  combustion.  In  actual  propellant  combustion, 
the  changes  of  the  chemical  equivalent  have  a  very  large  influence 
on  the  combustion  process.  We  still  await  test  proof.  The  reason 
for  the  production  of  super  burning  rate  combustion  and  plateau 
combustion  is  the  changes  in  the  course  of  ester  nitrate  decom¬ 
position.  It  is  produced  in  the  condensed  phase  and  is  therefore 
listed  in  the  condensed  phase  theory.  Yet  the  action  results 
occur  in  the  gas  phase’s  hissing  region  and  thus  can  also  be 
viewed  as  the  condensed  phase  -  gas  phase  theory. 

3.  The  Chelate  Theory 

Suh  and  Lenchitz  et  al  [6]  (1974)  used  detonation  tests  and 

at  the  same  time  determined  the  gas  product's  composition  method 
to  discover  the  50%  decrease  of  NO  in  the  combustion  products  of 
the  catalytic  propellant,  the  total  amount  of  non-changing 
CO 2  and  CO  and  the  C02/C0  ratio  which  actually  slightly  decreased. 
Simultaneously,  they  also  discovered  that  the  detonation  heat  of 
the  catalytic  propellant  rose  under  low  temperature  and  the 


combustion  surface  temperature  also  rose.  This  paper  considers 
that  this  was  the  result  of  the  increase  of  the  NO  reduction 
quantity  in  the  area  of  the  surface.  They  proposed  that  the 
lead  salt's  lead  atoms  and  the  oxygen  atoms  of  the  ester  nitrate 
form  a  chelate  type  compound  and  the  electrons  on  the  ester 
nitrate  shift  toward  the  divalent  lead.  This  causes  the  RO-NC>2 
chain  to  weaken  and  the  ester  nitrate  decomposition  quickens 
causing  the  burning  rate  to  increase. 

Early,  in  1966,  Lenchitz  et  al  [7,8]  researched  the  catalytic 
action  mechanism.  At  the  time,  the  ratio  of  the  measured  cata¬ 
lytic  propellant 's  combustion  product  C02/C0  increased.  This 
theory  explained  that  the  reason  why  the  covalent  compound  of 
tetraphenyl  lead  does  not  have  catalytic  action  is  because  it 
cannot  form  a  chelate.  Yet,  it  cannot  explain  lead  tetraethyl, 
lead  oxide  and  lead  powder  being  able  to  cause  double-base 
propellants  to  produce  plateau  combustion  [9],  There  is  no  direct 
evidence  that  the  chelate  accelerates  the  decomposition  of  the 
ester  nitrate  so  that  the  detonation  heat  increase,  when  surface 
temperature  rises  and  pressure  is  low,  can  also  result  in  even 
more  NO  being  reduced.  Reference  [10]  considers  that  the  reason 
the  double-base  propellant's  detonation  heat  drops  when  there  is 
low  pressure  is  the  expansion  of  the  pressurized  inert  gas. 

The  additional  nickel  powder  causes  the  detonation  heat  to 
greatly  increase  under  low  pressure  and  also  causes  the  luminous 
flame  to  be  close  to  the  combustion  surface.  The  nickel  powder, 
however,  does  not  have  a  plateau  reaction. 

4.  Spectral  Method  Research  on  the  Action  Position  and 
Mechanism  of  Catalysts 

Dr.  Mal'tsev  [11]  (1973)  showed  in  research  using  an  ultra¬ 

violet  transmitted  spectrum  that  the  catalytic  combustion  product 
CuH  in  the  area  of  the  double-base  propellant's  combustion 
surface  with  added  copper  salicylate  had  strong  absorption,  the 


dark  region  was  relatively  narrow,  the  N^O  dissolving  was  very 
fast  and  the  surface  temperature  rose.  The  copper  compound 
caused  the  CO  and  C02  to  appear  even  earlier  possibly  delaying 
the  appearance  of  1^0.  This  explains  the  accelerated  gas  phase 
reaction  and  thermal  feedback  toward  the  surface.  However,  the 
lead  salicylate,  the  temperature  along  the  flame's  tongues  of 
fire  direction  and  the  composition  had  no  noticeable  changes. 

The  lead  spectral  line  was  not  included.  This  also  explains  that 
the  lead  salacylate  greatly  accelerated  the  chemical  decomposi¬ 
tion  of  the  initial  condensed  phase  of  the  combustion  as  well  as 
the  action  of  the  solid  phase  transforming  into  a  gas  phase.  The 
lead  caused  the  H20  to  appear  even  earlier,  to  obtain  carbon 
bearing  substances  and  then  to  be  burned  off. 

Based  on  the  spectrum,  they  proved  that  the  lead  acted  in 
the  condensed  phase,  the  copper  acted  in  the  gas  phase  and  the 
two  acted  independently  in  different  combustion  regions.  Yet, 
the  later  mentioned  Androsov  et  al  considered  that  the  increased 
effect  on  the  lead  was  a  type  of  non-added  effect.  Baer  et  al 
[12]  of  the  United  States  also  confirmed  Mai 'tsev' s  results. 
Because  the  experiments  were  carried  out  in  actual  combustion 
systems,  they  were  more  believable.  Yet,  the  combustion  surface 
was  full  of  bumps  and  hollows  and  it  was  very  difficult  to 
differentiate  within  a  several  micrometer  to  several  ten  micro¬ 
meter  space  on  the  surface.  Therefore,  it  was  necessary  to 
raise  the  resolution  and  use  a  low  plateau  pressure  propellant 
in  order  to  be  able  to  further  show  the  chemical  catalytic 
nature. 

5.  Heat  Decomposition-Gas  Phase  Chromatographic  Research 
Results 

The  Androsov  Research  Group  used  heat  decomposition-gas 
phase  chromatography  as  the  basic  method  to  systematically  study 
the  action  mechanism  of  lead-copper  propellants.  On  the  one 
hand,  they  researched  the  characteristics  of  the  lead-copper 


composite  catalysis  such  as  the  influences  of  initial  tempera¬ 
ture  and  pressure  [13] ,  the  influence  of  the  propellant 
component  [14] ,  the  influences  of  the  lead-copper  ratio  [15] 
etc.  in  order  to  obtain  perceptual  knowledge.  On  the  other  hand, 
they  used  heat  decomposition  chromatography  to  find  chemical 
evidence  [16,17,18],  to  test  and  verify  the  obtained  perceptual 
knowledge  and  to  gain  rational  knowledge.  In  1978,  they 
published  "The  Mechanism  of  Lead-Copper  Composite  Catalysts  on 
Propellant  Combustion  Action." 

Their  research  showed  that  CuO  does  not  change  the  decom¬ 
position  speed  of  ester  nitrate  but  causes  the  amount  of  CO  in 
the  combustion  product  to  decrease,  the  amount  of  CO  to  increase 
and  the  CO  to  oxidize  into  CO^.  When  the  cellulose  nitrate  decom¬ 
poses,  it  is  necessary  to  have  the  participation  of  CuO  and  when 
the  glycerine  nitrate  decomposes,  it  is  necessary  to  have  the 
participation  of  N02»  The  action  of  the  lead  oxide  causes  the 
highest  quantity  of  heat  of  the  C02  product  in  the  decomposed 
product  to  increase.  It  also  causes  the  reaction  ability  ratio 
of  NO  and  N02  when  in  combustion  to  increase  the  quantity  of  the 
high  N20.  As  a  result,  super  burning  rate  combustion  is  produced. 
The  increase  in  strength  of  the  lead-copper  catalytic  action 
(exceeds  the  total  independent  action  of  the  two)  is  owing  to 
the  PbO  causing  the  C02  and  N20  to  increase  and  the  CO  to 
decrease  in  the  propellant's  decomposed  products.  Because  the 
decrease  of  CO  is  beneficial  to  the  catalytic  action  of  CuO, 
this  causes  even  more  CO  to  transform  into  C02- 

The  increase  of  glycerine  nitrate  in  the  propellant  causes 
the  action  of  the  CuO  to  drastically  drop  and  when  the  glycerine 
nitrate  exceeds  40%,  the  action  completely  disappears.  This  is 
because  the  CO  formed  when  the  glycerine  nitrate  decomposes  is 
much  greater  than  the  C02  and  the  surface  of  CuO  sustains  CO 
poisoning.  When  double-base  propellants  with  52%  of  diethyl 
glycol  ester  dinitrate  is  at  above  40  atm  pressure,  the  action 


of  CuO  is  stronger  than  the  action  of  PbO^  and  if  a  lead-copper 
composite  catalyst  is.  used,  the  catalytic  effect  is  greatly 
increased.  Yet,  when  in  70-100  atm  pressure,  mesa  combustion  is 
produced.  Diethyl  glycol  ester  dinitrate  decomposition  releases 
a  great  deal  of  ^0  and  produces  CO.  However,  the  N^O  causes 
the  burning  rate  to  increase  and  the  CO  causes  CuO  poisoning. 
This  indirectly  verified  the  action  mechanism  of  the  lead-copper 
composite  catalyst.  They  employed  common  procedures  to  rather 
fully  explain  the  chemical  qualities  of  the  production  of  the 
super  burning  rate  and  plateau  and  mesa  combustion.  Because  the 
heat  decomposition  and  combustion  processes  are  linked  as  well 
as  have  differences,  we  will  await  other  direct  test  verifica¬ 
tion. 


6.  Mass  Spectra  Research  Results 

Bauman  [20]  (1965)  and  Dauerman  [21,22]  (1968)  first  used 

mass  spectra  to  research  the  effects  of  catalysts  on  double¬ 
base  propellant  heat  decomposition.  Faber  [23]  (1976)  used  a 

specially  designed  four  differential  pump  mass  spectrometer  to 
study  the  products  formed  in  combustion  flames  under  constant 

pressure.  After  the  nozzle  dropped  in  pressure  many  times  the 
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speed  dropped  to  10  atm  pressure  and  this  molecular  beam 
entered  the  mass  spectrometer  for  analysis.  It  was  discovered 
that  when  lead-copper  oxyacid  salts  decompose,  metallic  lead 
and  copper  are  produced  but  the  non-catalytic  action  of 
manganese  salt  produces  MnO.  The  covalent  linked  tetraphenyl 
lead  vaporizes  but  does  not  decompose.  The  gaseous  lead  and 
copper  atoms  react  with  the  ester  nitrate  and  produce 
PbN02  and  CuN02.  These  compounds  are  stable  under  suitable  tem¬ 
peratures.  When  the  nitrite  enters  into  combustion  area,  a 
reaction  occurs  which  produces  hydrogen  oxides  and  oxides. 

This  is  consistent  with  thermal  dynamic  prediction.  Therefore, 
the  appearance  of  PbN02  and  CuN02  is  possibly  the  different¬ 
iating  standard  of  plateau  combustion. 


This  research  maintains  the  view  that  when  the  ester 
nitrate  decomposes,  the  CO-NC>2  separates  into  CO  and  N02  which 
does  not  verify  the  NO^  discovered  by  Bauman  and  Dauerman.  it 
does  not  coincide  with  Prekel's  Pb-PbO  cycle  oxidation  reduction 
reaction  but  is  identical  to  Suh's  hypothesis.  That  is,  the  N02 
uses  a  certain  method  and  the  lead  action  in  the  lead  salt  yet 
does  not  form  a  chelate.  It  does,  however,  form  PbN02>  The  gas 
state  Pb  and  N02  reaction  can  cause  the  ester  nitrate  to  release 
even  more  N02  and  thus  increase  the  burning  rate.  Mass  spectra 
is  a  powerful  tool  for  studying  condensed  phase  reactions.  If 
the  system  proposed  by  White  [25]  is  used  for  study  under  high 
pressure,  the  molecular  beam  will  instantly  be  sent  into  the 
mass  spectra  and  this  will  make  an  important  contribution  for 
the  plateau  catalytic  mechanism. 

7.  Results  of  Differential  Scanning  Calorimetric  Tests 

Eisenreich  et  al  [26,27]  (1976)  used  a  differential  scanning 

calorimeter  to  simultaneously  carry  out  heat  weightlessness  and 
heat  analysis  tests.  It  was  proven  that  the  first  step  of 
double-base  propellant  heat  decomposition  is  glycerine  nitrate 
vaporization  and  the  second  step  is  the  heat  decomposition  of 
the  cellulose  nitrate.  Under  low  pressure,  the  reaction  heat 
of  the  catalytic  propellant  must  be  higher  than  that  for  non¬ 
catalysis.  Its  law  is  identical  to  Lenchitz ' s  determinations 
using  the  detonation  heat  method  [8,9].  The  peak  temperature 
of  catalytic  propellant  decomposition  is  relatively  high  and 
active  energy  is  relatively  low.  Similar  work  done  domestically 
has  also  obtained  the  same  law. 

Eisenreich  et  al  [28]  (1978)  continued  to  research  the 
influence  of  copper  and  lead  compounds  on  cellulose  nitrate 
decomposition.  Beatriz  et  al  [29]  (1978)  then  used  a  different¬ 

ial  scanning  calorimeter  to  find  the  major  gas  phase  or 
multiphase  reaction  actions  of  double-base  propellant  trajectory 
modifiers  on  adjacent  combustion  surfaces.  The  heat  quantity 
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measured  by  the  differential  scanning  calorimeter  is  not  pure 
solid  reaction  heat  [30/31,32].  However,  Eisenreich  considered 
the  solid  phase  reaction  as  essential  and  the  catalyst  as 
primarily  acting  in  the  solid  phase.  Differential  scanning 
calorimetry  is  a  modern  heat  analysis  method  which  is  quite 
effective  for  researching  propellant  heat  decomposition.  Although 
its  heating  speed  differs  greatly  from  the  heating  speed  of  real 
propellants,  yet  its  results  still  have  reference  value.  The 
recently  developed  heat  analysis  method  which  uses  a  high  heating 
speed  (100-20,000  degrees/second)  is  employed  to  study  the 
ignition  dynamics  of  propellants  and  explosives  [33,34,35]. 

8.  Measuring  the  Detonation  Point  and  Finding  the 
Activation  Energy 

Bai  Mulan  et  al  (1978)  used  the  method  of  measuring  the 

detonation  point  to  find  the  reaction's  activation  energy  and 

index  forward  factor.  Under  200-250°C  temperature,  a  0.03  gram 

propellant  sample  is  used  to  measure  the  explosion  delay  period 

E 

under  6  temperatures.  When  making  the  InX  =  +  InC  curve,  we 

find  E  and  C  listed  in  the  table  as  follows  (in  the  formula: 

T  is  the  delay  period,  E  is  the  activation  energy,  R  is  the 
gas  constant  and  C  is  index  forward  factor) : 
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Key: 


1.  Calculation  results 

2.  Propellant 

3.  Activation  energy  E  (kilocalories/ 
gram  molecules) 

4.  Index  forward  factor  C  (1/second) 
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Key:  5.  Cellulose  nitrate  (I) 

6.  I  with  2.6%  lead  benzoate 

7.  I  with  2.6%  lead  tannate 

8.  A  certain  double-base  propellant  (II) 

9.  II  with  2.6%  lead  benzoate 

10.  II  with  2.6%  lead  tannate 

It  can  be  seen  from  the  data  listed  in  the  table  that  the 
catalyst  causes  the  double-base  propellant  and  cellulose  nitrate 
activation  energy  to  decline,  the  condensed  phase  reaction  to 
accelerate  and  the  burning  rate  to  increase. 

The  activation  energy  measured  by  this  method  is  naturally 
the  explosion  reduction  period’s  activation  energy  which  is 
different  from  the  concept  of  explosion  reaction  or  combustion 
reaction  activation  energy.  If  there  is  identical  control  of 
the  reaction,  the  activation  energy  will  also  be  the  same. 

Zenin  [36]  proposed  a  model  for  the  double-base  propellant's 
condensed  phase  reaction  layer.  He  proved  that  the  activation 
energy  of  ester  nitrate  decomposition  is  about  48-50  kilocalories  / 
gram  molecule  and  the  combustion  activation  energy  is  19  kilo¬ 
calories/gram  molecule  activation  energy  of  the  N02  oxidized 
hydrocarbon.  Therefore,  the  control  reaction  of  double-base 
propellant  combustion  is  possibly  the  oxidation  reaction  of  NC>2. 
The  activation  energy  measured  by  the  detonation  point  method  is 
basically  the  same  as  the  combustion  activation  energy.  It  is 
very  possibly  the  same  as  the  control  reaction,  that  is,  the 
oxidation  reaction  of  NC^.  The  reaction  when  measuring  the 
detonation  point  is  not  a  pure  solid  phase  reaction  but  should 
include  the  gas  phase  reaction  of  the  solid  phase  surface.  This 
can  use  simple  test  proof. 

9.  Another  Result  in  Measuring  the  Surface  Temperature 

Denisyuk  .  et  al  [18]  (1971)  used  a  thermal  couple  to 
measure  the  temperature  distribution  curve  of  catalytic  and  non- 
catalytic  H  type  explosives  under  20-90  atm  pressure.  Moreover, 
they  derived  data  on  the  surface  temperature,  thickness  of  the 


stack  gas  area,  the  highest  flame  temperature,  the  distance 
from  the  surface  to  the  highest  flame  temperature,  the  heat 
released  by  the  condensed  phase,  the  heat  released  by  the 
stack  gas  area  and  the  temperature  gradient  in  the  vicinity  of 
the  surface.  They  discovered  that  the  catalyst  did  not  change 
the  distance  from  the  surface  to  the  highest  flame  so  that  the 
gas  phase  heat  transfer  decreased  4-5%  and  the  stack  gas  area's 
heat  loss  also  decreased.  Therefore,  the  catalyst  must  act  in 
the  condensed  phase  to  cause  the  burning  rate  to  increase.  Under 
high  pressure,  the  gas  product's  solubility  rises  in  the  con¬ 
densed  phase  which  causes  catalytic  poisoning  or  the  catalyst  to 
become  an  inhibitor.  As  a  result,  the  burning  rate  drops 
producing  plateau  and  mesa  combustion. 

I  used  the  thermal  couple  to  measure  the  temperature  so  as 
to  prove  that  catalysis  occurs  in  the  condensed  phase,  yet  I  do 
not  know  the  concrete  chemico-physical  process.  Under  high 
pressure,  it  is  difficult  to  measure  the  surface  temperature 
correctly.  The  fact  that  the  catalyst  does  not  change  the 
distance  from  the  surface  to  the  highest  flame  is  opposite  of 
Kuboto's  measurement  results  mentioned  in  section  2. 

III.  The  Gas  Phase  Theory 

1.  The  Lead-Lead  Oxide  Cycle  Theory 

Steiberger  [1]  (1952)  added  organic  lead  salt  into  a  propel¬ 

lant  after  decomposing  it  in  a  furnace  and  discovered  that  the 

o 

catalytic  effects  were  better.  When  in  350  C,  the  heat  decom¬ 
posed  product  was  lead  or  lead  oxide.  They  should  be  the  cause 
of  the  propellant's  plateau  combustion.  The  organic  part  of 
the  lead  salt  only  acted  to  produce  fine  metallic  lead  or  lead 
oxide.  For  this  reason,  he  proposed  the  gas  phase  catalytic 
mechanism: 


Pb+NO,(NO) — ►PbO 
Pb0+C0(CH,0,  H»,  Cm) — *Pb 
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Prekel  137]  (1961)  used  a  sealed  exploder  to  measure  the 
burning  rate.  He  observed  a  broken  piece  of  carbon  fall  down 
in  the  plateau  pressure  range  and  there  was  a  melted  lead  ball. 
This  phenomenon  did  not  occur  above  the  plateau  pressure.  He 
I  considered  that  the  above  mentioned  reaction  occured  in  the  gas 
phase  area  or  foam  area.  Preckel  [38]  (1965)  also  explained, 

based  on  the  Pb-PbO  cycle  theory,  that  inorganic  lead  compounds 
and  aliphatic  lead  compounds  can  only  cause  a  low  heat  value 
catalyst  to  produce  plateau  combustion.  However,  aromatic  com¬ 
pounds  can  cause  a  high  heat  value  catalyst  to  produce  plateau 
combustion. 

This  theory  is  quite  rough.  The  lead- lead  oxide  cycle 
reaction  is  indirectly  infered  and  is  not  directly  evidenced. 
Preckel  indicated  the  area  of  catalytic  action  and  discovered 
that  the  broken  piece  of  carbon  was  related  to  plateau  combus¬ 
tion.  This  opened  avenues  for  the  lead  and  lead  containing 
substances  theory. 

2.  The  Free  Base  Theory 

Sinha  et  al  [39]  (1968)  consider  that  the  burning  rate  of 

double-base  propellants  is  determined  by  the  extinction  speed 
and  the  extinction  speed  sustains  the  control  of  the  hissing 
region  on  the  surface  heat  transfer  but  is  not  influenced  by 
the  flame  area.  They  adopted  the  test  results  of  Paneth  et  al 

regarding  the  ability  of  lead  to  react  with  a  free  base  at 

400°C  to  produce  lead  alkyl  and  the  ability  of  lead  alkyl  to 
decompose  at  600°C  into  lead  and  a  free  base.  The  lead  and 

ester  nitrate  decomposed  free  base  combine  in  the  foam  area  to 

form  lead  alkyl  and  in  the  hissing  area  the  lead  alkyl  decom¬ 
poses  into  a  free  base  and  lead.  The  free  base  produces  a  chain 
reaction  which  causes  NO  reduction,  emits  a  large  amount  of 
heat,  accelerates  the  transfer  heat  from  the  hissing  area 
towards  the  foam  area,  causes  the  burning  rate  to  rise  and 
produces  super  burning  rate  combustion.  However,  for  the  non- 
catalytic  propellant,  the  NO  does  not  participate  in  the 


reaction  of  the  hissing  area.  When  the  pressure  rises,  the 
free  base  concentration  in  the  hissing  area  rises  and  the  free 
base  automatically  stops  the  collision  increase  which  decreased 
the  reaction  between  the  free  base  and  NO  and  lowers  the 
transfer  heat  from  the  hissing  area  to  the  foam  area.  On  the 
other  hand,  when  the  pressure  rises,  the  heat  induction  coef¬ 
ficient  of  the  gas  increases  and  the  heat  transfer  increases 
from  the  hissing  area  towards  the  foam  area.  In  a  certain 
pressure  range,  these  two  opposite  effects  are  equal  and  this 
produces  plateau  combustion. 

The  free  base  theory  of  cellulose  nitrate  decomposition  was 
proposed  early  by  Wolfram  et  al  [40].  Morrow  [41]  (1976)  used 

a  heat  platform  polarizing  microscope  to  observe  the  lique¬ 
faction  decomposition  process  of  the  pure  cellulose  nitrate 
and  modified  cellulose  nitrate  (adding  a  free  base  inhibitor, 
accelerating  agent,  burning  rate  catalyst  and  additive)  when 
heated  to  225°C.  Test  results  showed  the  occurence  of  an  import¬ 
ant  free  base  process.  This  theory  did  not  study  the  free  base 
process  in  actual  propellant  combustion  and  thus  cannot  explain 
the  mesa  effect.  Moreover,  a  large  number  of  non-lead  metals 
possibly  react  with  the  free  base  and  produce  alkyl  metal  free 
base  yet  there  is  no  plateau  combustion. 

3.  The  Carbon  Bearing  Substance  -  Hot  Light  Sphere  Theory 

Eisenreich  et  al  [42]  (1978)  used  micro-high  speed  photography 
and  a  scanning  electron  microscope  to  observe  the  color  of  the 
solid  phase  area  40  micrometers  below  the  combustion  flame-out 
surface  change  from  a  certain  transparency  to  a  brown  or  brownish 
red  color.  They  discovered  that  the  largest  diameter  was  a  hollow 
of  50  micrometers  and  carbon  covered  the  entire  flame-out  surface. 
The  carbon  particles  had  20-50  micrometer  coarse  branches  and  the 
largest  diameter  was  a  100  micrometer  flat  plate.  The  carbon 
particles  were  porous  and  the  diameter  of  the  holes  were  about 
5  micrometers.  When  the  carbon  particles  were  above  5 


atmospheric  pressure  there  began  to  appear  a  red  light  which 

seemed  to  be  the  result  of  being  oxidized  by  the  decomposed 

products  flowing  by.  At  the  same  time,  the  carbon  particles 

changed  little,  they  flew  from  the  surface  and  their  number 

decreased.  The  diameters  of  the  carbon  particles  in  the  plateau 

pressure  range  were  about  20-50  micrometers  and  density  was  300 
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particles/centimeter  .  When  the  pressure  exceeded  5  atmospheric 
pressure,  metallic  light  spheres  appear  on  the  carbon  particles 
or  combustion  surface.  The  diameter  is  about  100  micrometers, 
the  light  spheres  which  flew  from  the  surface  have  diameters  of 
about  20-50  micrometers  and  their  speed  is  15-50  centimeters/ 
second.  They  consider  that  the  trajectory  modifier  accelerates 
the  decomposition  of  cellulose  nitrate  below  the  combustion 
surface  when  under  low  pressure  and  catalyzes  the  oxidation  of 
carbon  bearing  substances  on  the  combustion  surface.  As  a 
result,  a  great  deal  of  energy  is  released  and  super  burning 
rate  combustion  is  produced.  Because  the  carbon  bearing  and 
metal  bearing  light  spheres  leave  the  combustion  surface  and 
cause  the  energy  conducted  to  the  surface  to  diminish,  plateau 
and  mesa  combustion  are  produced. 

Early,  in  1956,  when  Brown  et  al  [1]  used  micro-high  speed 
photography  to  observe  the  structure  of  the  plateau  propellant' 
combustion  surface  they  discovered  that  under  plateau  pressure 
there  are  30-60  micrometer  diameter  lead  spheres.  In  the  last 
several  years,  a  good  deal  of  domestic  research  work  has  con¬ 
firmed  the  carbon  bearing  substance  -  hot  light  sphere  theory. 
It  can  directly  describe  the  forms,  sizes,  number  and  movements 
of  the  carbon  bearing  substances  and  metallic  light  spheres 
related  to  the  production  of  plateau  combustion.  Although  its 
chemical  properties  have  not  yet  been  precisely  proven,  yet  its 
possible  mechanism  has  been  suggested.  Reactions  on  the  carbon 
particles  have  confirmed  that  the  carbon-lead  catalysis 
mentioned  below  produces  a  mechanism  of  super  burning  rate 
combustion. 


4.  The  Lead-Carbon  Catalysis  Theory 

Hewken  et  al  [43]  (1971)  published  a  summary  of  the  work  of 

their  research  team.  They  used  a  flame  probe  to  measure  the 
composition  of  each  combustion  area,  used  a  heat  couple  and  the 
infrared  heat  measuring  method  to  measure  the  temperature  bf 
the  combustion  surface  and  used  high  speed  photography  to  study 
the  structure  of  the  combustion  waves.  They  discovered  that  under 
low  pressure,  the  introduction  of  a  lead  compound  caused  the 
temperature  near  the  surface  to  rise;  analysis  of  the  catalytic 
propellant's  dark  area  showed  that  the  ^  in  the  combustion 
products  greatly  increased  and  the  NO  decreased;  there  was  a 
great  amount  of  heat  released  from  the  NO  reduction  reaction; 
the  action  of  the  lead  compound  caused  the  combustion  surface  to 
produce  a  large  amount  of  carbon.  There  was  a  great  amount  of 
carbon  when  in  super  burning  rate  combustion  and  there  was  carbon 
produced  on  the  combustion  surface  within  the  whole  pressure 
range  yet  the  amount  of  carbon  was  much  less  under  low  tempera¬ 
ture  than  the  catalytic  propellant.  They  considered  that  the 
action  of  the  lead  compound  not  only  increased  the  production  of 
carbon  but  also  caused  carbon  activation.  The  activated  carbon 
is  the  solid  phase  catalyst  of  the  NO  released  heat  reduction 
reaction  and  when  there  is  a  non-catalyst,  the  speed  of  the  NO 
in  the  dark  area  reduction  reaction  is  very  slow.  As  a  result, 
the  heat  of  the  combustion  surface  increases  and  the  burning 
rate  also  increases.  Yet,  the  lead  oxidizes  the  NO  and  carbon 
reaction  causing  the  carbon  to  disappear.  When  the  pressure 
rises,  the  carbon  disappears  more  easily,  the  NO  reduction 
reaction  slows  down  and  the  propellant's  buring  rate  declines 
which  produces  mesa  or  plateau  combustion. 

Hewkin  et  al  [44]  (1974)  continued  research  on  the  connection 
between  combustion  surface  phenomena  and  the  catalyst  in  order  to 
prove  the  above  mentioned  theory.  The  observed  copper  compound 
had  a  slower  carbon  eliminating  speed  than  the  lead  compound  and 
therefore  the  catalysis  maintained  high  pressure.  When  the  lead 


and  copper  acted  simultaneously,  the  elimination  of  carbon  was 
determined  by  the  lead.  By  adding  an  oxygen  rich  substance  of 
ammonium  perchlorate  or  glycerine  nitrate  the  carbon  on  the 
surface  greatly  decreased  and  the  plateau  effect  was  destroyed. 
The  joining  of  the  cellulose  nitrate  accelerated  the  production 
of  carbon  which  caused  the  plateau  effect  to  expand  to  even 
higher  pressure.  Hexogen  does  not  change  the  flame  structure. 
There  is  still  a  lead  carbon  covering  yet  the  plateau  effect  is 
destroyed.  This  is  possibly  because  the  hexogen's  decomposed 
product  NO ^  causes  catalyst  poisoning.  It  explains  that  although 
the  formation  of  carbon  is  important  the  reaction  is  also  very 
fast  [45]  at  200°C  and  the  simultaneous  reaction  of  copper  and 
lead  is  easier  to  carry  out  [46].  Therefore,  the  catalysis  of 
the  lead  and  copper  themselves  should  be  consiered  the  same  as 
the  carbon. 

This  theory  can  explain  many  phenomena.  However,  although 
the  reaction  of  lead  oxidizing  NO  and  carbon  consumes  carbon 
the  amount  of  released  heat  is  very  great  and  thus  cannot  offer 
a  tenable  argument.  It  appears  that  the  use  of  Eisenreich's 
research  can  explain  plateau  combustion.  That  is,  the  carbon 
bearing  and  metal  bearing  light  spheres  leaving  the  surface 
cause  the  energy  transmitted  to  the  surface  to  diminish  and 
produce  plateau  combustion. 

5.  The  7T  Chain  Complex  Theory 

Fifer  et  al  [47]  (1975)  reported  that  when  a  laser  ignited 

a  thin  film  sample  of  cellulose  nitrate  and  double-base  pro¬ 
pellant  and  it  was  burned  under  1,  20  and  40  atm  pressure,  an 
infrared  spectrum  was  used  to  measure  the  final  product.  Under 
super  burning  rate  combustion  pressure,  a  lead  salt  bearing 
propellant  sample  obtained  an  increase  in  the  CC^/CO  ratio  in 
the  combustion  product.  This  caused  the  surface  and  dark  area 
temperatures  to  increase  and  the  burning  rate  to  increase. 

There  was  no  notice. -le  effect  of  the  water  and  NO.  permeating 


into  the  _est  sample  on  the  combustion  product.  However,  a 
small  amount  of  02  catalyzed  a  reaction  in  which  the  CO 
oxidized  into  C02  and  the  NO  reduced  to  N2*  They  used  Salooja's 
[48,49]  (1967)  7C  chain  complex  theory.  That  is,  the  PbO  inhibits 

the  hydrocarbon  combustion  yet  forms  a  Jl  chain  surface  complex 
with  an  oxygen  bearing  hydrocarbon  derivative.  This  greatly 
accelerates  the  combustion  of  the  oxygen  bearing  derivative. 

The  special  features  of  the  catalytic  reaction  are  that  the 
amount  of  C02  produced  increases,  sometimes  the  amount  of  H 
produced  also  increases,  the  PbO  is  reduced  to  Pb  and  the  Pb  is 
then  oxidized  into  PbO.  The  PbO  uses  the  same  type  of  initial 
stage  and  quadratic  solid  phase  decomposition  reaction  for  the 
catalysis  of  the  ester  nitrate. 

•  In  his  experiments,  Fifer  discovered  that  there  was  too 
little  carbon  which  was  not  sufficient  to  cause  a  noticeable 
change  in  the  N02/CH20  ratio.  As  a  result,  he  did  not 
support  Kubota's  chemical  equivalent  shift  theory.  Their 
measured  NO  reduction  quantity  was  very  small  and  thus  did  not 
coincide  with  the  large  amount  of  NO  being  catalyzed  by  the 
reduced  carbon- lead  and  the  chelate  theory.  However,  test 
results  were  based  on  measuring  laser  ignited  thin  film  com¬ 
bustion  products  which  even  though  differing  from  real 
propellant  combustion  were  still  worth  studying. 

6.  Research  Results  of  the  Hot  Platform  Microscope 

Morrow  (1973-1974)  used  a  hot  platform  microscope  in  research¬ 
ing  thermal  high  molecular  degradation  as  well  as  for  a  great  deal 
of  work  in  the  areas  of  the  condensed  phase  heat  properties  of 
cellulose  nitrate  with  different  degrees  of  nitrate  and  combustion 
propellants  [50,51,52].  He  specially  researched  the  influence  of 
catalysts  on  cellulose  nitrate  thin  film  combustion  [53,54],  A 
180-40  micrometer  thin  film  sample  in  a  capillary  melting  point 
analyzer  with  an  outer  diameter  of  1.5  millimeters  was  heated  to 
250°C,  a  polarizing  microscope  was  used  to  observe  the  phenomena 
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during  the  heating  process  and  the  heat  loss  in  weight  was 
measured.  Four  sets  of  tests  were  carried  out  on  inorganic  and 
organic  lead  salicylate  using  1  to  34  atm  pressure  and  it  was 
discovered  that  there  was  no  great  difference  between  the 
appearance,  beginning  bubble  and  ending  bubble  temperatures  of 
the  flower  shaped  substance.  However,  the  amount  of  weight  lost 
was  opposite  of  the  prediction  and  the  loss  of  weight  of  the 
added  catalyst  was  actually  quite  small.  He  considered  that  the 
catalyst  does  not  act  upon  the  solid  phase  reaction  but  only 
acts  on  the  gas  phase  above  the  surface.  This  is  identical  to 
Hewkin 's  view.  It  appears  that  this  type  of  method  can  only  aid 
the  procedure  but  cannot  basically  resolve  the  problem. 

7.  Other  Evidence  for  the  Gas  Phase  Theory 

Lengelle  et  al  [55]  (1979)  used  the  concept  of  a  composite 
propellant's  combustion  mechanism  to  study  the  mechanism  of  a 
double-base  propellant's  plateau  combustion.  They  synthesized 
their  own  test  data  with  that  of  past  references  and  proved 
that  the  surface  temperature  and  combustion  assumes  a  linear 
relationship  and  that  the  catalyst  does  not  change  this  rela¬ 
tionship.  The  condensed  phase  heat  emission  quantity  calculated 
from  the  surface  temperature  greatly  increases  with  the  increase 
of  the  burning  rate  yet  is  not  related  to  the  catalyst.  There¬ 
fore,  the  propellant's  decomposition  dynamics  are  not  influenced 
by  the  catalyst,  the  decomposition's  manifested  activation 
energy  is  abouu  40  kilocalories/gram  molecule  and  is  possibly 
correlated  with  the  C0-N02  chain  breaking.  When  a  porous  metallic 
pyrolytic  instrument  was  used  to  study  the  decomposition  of 
propellants  in  a  vacuum  and  its  produced  gas  was  determined  by 
a  mass  spectrometer,  it  was  discovered  that  the  influence  of  the 
catalyst  on  the  gas  composition  was  not  great.  Further,  as 
regards  decomposition  dynamics,  the  copper  and  lead  salt  has  still 
not  changed  in  the  condensed  phase.  A  conclusion  is  made  from  this 
that  although  the  propellant  releases  a  large  amount  of  heat  in 
the  condensed  phase  the  catalyst  does  not  act  upon  the  condensed 


phase  process.  The  active  products  PbO  and  CuO  of  the  catalyst 
only  act  in  the  initial  flame  and  therefore  the  precision 
properties  and  stability  of  the  lead  salt  are  not  major  factors. 
They  used  the  reaction  of  PbO  and  formaldehyde  in  a  simulation 
reactor  (residence  time  was  about 0.1  seconds,  the  temperature 
reached  350°C)  to  prove  that  this  reaction  very  possibly  pro¬ 
duces  carbon  on  the  propellant's  surface.  They  also  discovered 
that  the  end  temperature  of  the  initial  flame  rises  with  the 
pressure  and  when  there  is  super  burning  rate  combustion  the 
corresponding  initial  flame  end  temperature  is  even  higher. 
Therefore,  the  action  position  of  the  catalyst  is  in  the  gas 
phase's  initial  flame  area.  When  there  is  lead  and  copper  salt, 
the  PbO  and  formaldehyde  reaction  in  the  initial  fame  produces 
carbon  dregs  causing  even  more  NO  to  participate  in  the  reaction 
and  thus  releases  more  heat  and  produces  super  burning  rate 
combustion. 

A  porous  plate  pyrolytic  instrument  was  used  to  study  the 
decomposed  products  which  disappeared  in  the  condensed  phase 
when  there  was  no  gas  phase  reaction.  When  research  was  already 
done  on  propellants  with  no  catalyst,  then  when  there  is  a 
catalyst  research  should  be  continued.  Evidence  of  gas  phase 
catalysis  is  obtained  indirectly  and  we  still  await  research 
for  direct  evidence.  The  method  of  researching  composite 
propellants  is  very  useful  and  should  be  extended  to  many  other 
items . 

8.  The  Theory  of  the  Reaction  With  No  Relation  Between 
the  Hissing  Area  and  Pressure 

In  order  to  obtain  more  accurate  surface  temperature  and 
clearer  photographs  of  each  area  of  the  flame,  Kubota  [56] 

(1979)  used  a  specially  developed  formulation  with  low  plateau 
pressure,  a  wide  range  and  a  low  burning  rate.  This  plateau 
propellant  had  very  strong  super  burning  rate  combustion 
properties  under  16  atm  pressure  and  under  16-36  atm  pressure 
the  plateau  burning  rate  was  0.24  centimeters/second.  The 
corresponding  non-catalytic  propellant's  pressure  index  was  0.85. 


53 


Research  found  that  the  dark  areas  total  reaction  progression 
was  2.6  which  is  unrelated  to  the  existence  of  the  catalyst. 

This  shows  that  the  catalyst  does  not  change  the  process  of  the 
dark  area.  The  burning  rate  does  not  change  in  the  plateau  area 
but  increases  with  the  pressure.  The  catalytic  propellant's 
luminous  flame  area  approaches  the  combustion  surface  faster 
than  the  non-catalytic  propellant's  luminous  flame  area.  The 
temperature  gradient  of  the  non-catalytic  propellant's  hissing 
area  uniformly  rises  with  the  pressure  but  the  temperature 
gradient  of  the  catalytic  propellant's  hissing  area  rises  in 
the  super  burning  rate  region  and  does  not  change  in  the  plateau 
region.  The  hissing  area  temperatures  of  the  two  types  of 
propellants  change  with  the  pressure  and  are  the  same  as  the 
burning  rate  law  of  changing  with  the  pressure.  The  temperature 
gradient  of  the  catalytic  propellant's  hissing  area  remains 
invariant  in  the  plateau  region,  that  is,  the  reaction  speed  does 
not  change.  As  a  result,  the  feedback  of  the  heat  conducted  from 
the  hissing  area  to  the  combustion  surface  remains  invariant  in 
the  entire  pldteau  region.  Nevertheless,  the  catalytic  activity 
in  the  plateau  area's  hissing  region  is  changing  and  moreover  it 
rapidly  decreases  with  the  pressure  increase.  This  is  a  result 
of  the  related  inhibited  reaction  of  the  hissing  area  and  the 
lead.  The  total  reaction  progression  of  Kubota  [1]  and  the 
results  of  Suh  [6]  and  Lengelll  [55]  but  is  identical  to 
Oenisyuk's  results  [18].  However,  based  on  the  observation  of 
Denisyuk  ,  the  catalyst  does  not  change  the  distance  from  the 
surface  to  the  highest  flame.  This  is  opposite  to  Kubota's 
observations.  Thus,  Denisyuk  obtained  the  condensed  phase 
catalysis  theory  and  Kubota  obtained  the  gas  phase's  hissing 
area  catalysis  theory.  Because  of  this,  use  of  surface  tempera¬ 
ture  and  combustion  area  photographs  to  study  the  plateau 
catalyst's  action  mechanism  still  requires  continued  work.  We 
can  thus  hope  for  uniformly  recognized  conclusions. 


IV.  Several  Views  on  Future  Research  Directions 


1.  Research  Ohjects 

We  should  turn  from  simulated  real  combustion  systems  to 
real  combustion  systems.  Here,  real  combustion  systems  point  to 
real  propellant  formulations  for  combustion  under  high  pressure 
but  does  not  refer  to  combustion  in  real  engines.  To  date,  aside 
from  high  speed  photography,  surface  temperature  tests,  discon¬ 
tinuing  tests  etc.  on  physical  research  methods  and  flame  infrared 
emission  spectrum,  research  has  differed  greatly  with  real  pro¬ 
pellant  combustion  conditions.  In  this  way,  only  indirect,  refer¬ 
ence  type  data  can  be  obtained  so  that  definite  conclusions 
cannot  be  made.  Half  the  results  are  obtained  with  twice  the 
effort  and  there  is  a  good  deal  of  waste.  Therefore,  in  the 
future,  tests  which  depart  from  real  combustion  conditions  must 
be  avoided.  Naturally,  relevant  simulated  tests  are  still 
indispensable  supplementary  measures. 

2 .  Test  Methods 

Using  the  methods  of  high  speed  microphotography,  heat  couple 
heat  measuring  and  observations  of  the  flame-out  surface, 
Eisenreich  and  Kubota  have  reached  high  levels  in  their  research. 
They  obtained  a  large  number  of  active  physical  phenomena  of  the 
plateau  effect  and  inferred  the  chemical  properties  of  the 
plateau  effect  from  the  physical  phenomena.  For  example,  they 
inferred  the  possible  chemical  reactions  from  observation  research 
on  the  carbon  bearing  substances,  lead  bearing  substances  and 
light  spheres  on  the  combustion  surface.  Further,  they  obtained 
the  pressure  index  of  each  area  and  the  chemical  reaction 
progression  and  catalytic  activity  coefficient  by  measuring  the 
temperature  distribution  and  thickness  of  each  area  of  combustion 
and  carried  out  basic  judgements  for  determining  the  plateau 
catalyst's  action  position  and  the  chemical  reaction  types.  This 
technique  should  be  further  improved.  For  example,  raise  the 
test  accuracy,  develop  the  measuring  of  the  flame's  gas  phase 
temperature  distribution  and  use  color  high  speed  photography. 
Flame-out  surface  research  should  also  be  changed  from  physical 


phenomena  to  chemical  components.  There  are  some  units  in 
China  which  are  expanding  the  use  of  electronic  probes  to 
research  a  propellant’s  flame-out  surface. 

Mass  spectra  are  an  important  measure  for  studying  condensed 
phase  dynamics  which  should  be  given  serious  attention.  The 
joint  use  of  the  molecular  beam  sampling  technique  and  time 
jumping  mass  spectra  can  possibly  be  employed  to  study  the 
instant  products  of  flames  under  real  combustion  conditions.  This 
is  extremely  useful  for  solving  the  plateau  combustion  mechan¬ 
ism.  Naturally,  there  are  still  many  difficult  problems  which 
must  be  overcome  [25] . 

Laser  holography  and  laser  Raman  spectrum  are  combustion 
examination  techniques  which  await  development  [57,58].  Their 
advantage  lies  in  examining  real  combustion  systems  under  non¬ 
interference  combustion  conditions.  Laser  holographs  can  also 
measure  flame  temperature  distribution  and  observe  the  combus¬ 
tion  surface 's  physical  structure,  th^  catalysts  etc.  above  the 
surface  and  the  size,  distribution  and  movement  speed  in  the 
flame  when  existing  in  a  flame.  The  laser  Raman  spectrum  can 
measure  the  temperature  of  each  point  in  the  flame  and  the 
physical  types  and  concentrations. Naturally ,  it  is  a  superior 
examining  technique  for  combustion  chemical  reactions.  Among 
them,  because  the  coherent  anti-Stokes  Raman  scattering  is  very 
strong,  can  be  locally  analyzed  and  does  not  have  fluorescent 
interference,  it  can  be  considered  first  for  development. 

3.  Test  Designs 

Even  if  scientific  and  ingenious  design  tests  use  common 
test  methods  they  can  still  obtain  new  and  important  conclusions. 
For  example,  Kubota  designed  a  propellant  for  plateau  combustion 
under  16-30  atm  pressure  wherein  photography  and  temperature 
measuring  accuracy  were  greatly  increased  and  new  reliable 
conclusions  were  obtained.  Further,  in  order  to  verify  the  effect 
of  the  CH_0  and  NO,  reaction  in  the  propellant,  use  of  a 


carefully  designed  formation  can  be  envisioned  for  research. 

There  are  still  bright  prospects  in  the  area  of  test  designs. 

4.  Other  Problems  Related  to  Research  and  Plateau 
Combustion 

For  example:  the  effect  of  carbon  on  plateau  and  most  double¬ 
base  propellant  combustion;  the  mechanism  of  ultraviolet  light 
absorbers  etc.  raising  the  burning  rate  in  the  plateau  area;  the 
auxiliary  catalytic  effect  of  copper  on  lead;  the  effect  of 
inertial  additives  in  double-base  propellants;  non-lead  type 
plateau  catalysts  etc.  These  problems  are  all  related  to  the 
plateau  combustion  mechanism. 

To  sum  up,  although  the  plateau  combustion  mechanism  of 
double-base  propellants  has  been  researched  by  many  scientists 
throughout  the  world  for  over  20  years,  the  various  theories 
and  explanations  proposed  still  await  truly  substantial  knowledge 
of  this  problem.  Further,  a  large  amount  of  research  work  still 
awaits  us. 
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